Ancilla models for quantum operations: 
For what unitaries does the anciUa state have to be physical? 

Zhang Jiang/ Marco Piani,^ and Carlton M. Caves^'^ 

^ Center for Quantum Information and Control, University of New Mexico, 

MSC07-4220, Albuquerque, New Mexico 87131-0001, USA 

'^Institute for Quantum Computing and Department of Physics and Astronomy, 

University of Waterloo, Waterloo, ON N2L 3C1, Canada 

^Centre for Engineered Quantum Systems, School of Mathematics and Physics, 

The University of Queensland, St. Lucia, Brisbane 4072, Australia 

(Dated: March 2012) 

Cn , Any evolution described by a completely positive trace-preserving linear map can be imagined 

^~^ • as arising hom the interaction of the evolving system with an initially uncorrelated ancilla. The 

>— ' ' interaction is given by a joint unitary operator, acting on the system and the ancilla. Here we 

^^ , study the properties such a unitary operator must have in order to force the choice of a physical — 

^ . that is, positive — state for the ancilla if the end result is to be a physical — that is, completely 

Q ' positive — evolution of the system. 

r^ ■ I. INTRODUCTION 

(N; 

I— 1 1 Consider a primary system A, with Hilbert space Ha of dimension cIa, which is subjected to a superoperator £, 

i-G . and let B be an ancillary system, with Hilbert space Hb of dimension ds ■ It is a commonplace that if the action of 
^~-*' £ can be written as 

S : £{pa) - trs {UpA (E> ctbU^) , (1.1) 

C^ . 

^ ■ where [/ is a joint unitary operator acting on A and B and ctb is a physical state of B (a unit-trace, positive operator, 

^^] i.e., a density operator), then £" is a trace-preserving, completely positive (CP) map and is called a trace-preserving 

quantum operation [ll, I4I. Moreover, any quantum operation can be written in the form (|l.ll) for some joint unitary U 

^<\ ' and some ancilla state as- Quantum operations are thus the superoperator maps that can be realized physically [JS]. 

^ : The form (jl.ip for a quantum operation is sometimes called a Stinespring extension [J| , and we call it here an ancilla 

*/~J model. 

^^ , In this paper we are interested in a different question. Suppose we know that £, written as in Eq. (jl.ip with U a 

unitary operator, is a trace-preserving quantum operation. Can we conclude that as is a density operator? Generally 
not, as the case of a unitary transformation Ua of A makes immediately clear. Then we have U = Ua ® Ib and 






^l-j ■ £{pa) = ^t:b{<^b)UaPaUa^ so ctb can be any unit-trace operator. A natural question asks for the conditions on U 

~^ such that £ being a trace-preserving quantum operation implies that (Jb is a density operator. Here we formulate and 

I • prove a theorem that answers this question. 

• • I In Sec. ini we spell out the problem clearly and introduce mathematical concepts that are useful in addressing the 

. !^ ■ problem and proving our main result. Section IIIII states and proves our theorem, and Sec. IIVI explores a number of 

^^ ■ examples that illustrate features of the main result. Section FVl demonstrates a curious connection to the notion of 

^ , indirect state tomography, i.e., tomography on the system that can determine the state of the ancilla, and Sec. IVII 

_Cy_' concludes. 

II. MATHEMATICAL BACKGROUND 

Before stating our theorem, we introduce some mathematical notation and concepts. Consider a system Q with 
Hilbert space T-Lq of dimension dg. We will consistently refer to arbitrary elements of T-Lq as vectors and to vectors 
normalized to unity as state vectors. We refer to density operators, i.e., unit-trace, positive operators, as states; pure 
states are rank-one density operators, i.e., one-dimensional projectors corresponding to state vectors. We sometimes 
use a bra-ket notation for operators on T-Lq, i.e., O = \0) and O^ = (0|, and we denote the Hilbert-Schmidt operator 
inner product by {N\0) = iT:{N'^0). We use a tilde to distinguish operators that are normalized to unity, i.e., 
(0|0) = tr(Oto) ^ 1. 

A superoperator is a linear map on operators. Any superoperator Tq can be written as 

^Q=^^al3TaQTl^^Tal3\Ta){Tfs\ , (2.1) 



where the operators fa make up an orthonormal operator basis. The ordinary action of J-q on an operator O 
corresponds to replacing the by O: 

•^Q(O)-^^a0T„Of^. (2.2) 

The left-right action corresponds to using the bra-ket form of J-q to operate to the left or right on operators. A 
superoperator Tq has a Kraus decomposition [5| if it can be written as 

•^Q = E^"®^"=El^")(^"l- (2.3) 

a a 

The operators Ka are called Kraus operators. Having a Kraus decomposition is equivalent to J-q being positive in its 
left-right action. 

The standard way of defining complete positivity of a superoperator invokes a reference system R. A superoperator 
J-Q is completely positive if Tr ® Tq maps positive operators to positive operators for all reference systems R, 
where Z^ = IrQ Ir is the identity superoperator on R. It is clear that if J^q is left-right positive, i.e., has a Kraus 
decomposition, Tr®J-q maps positive operators to positive operators, thus making J-q completely positive. Moreover, 
if R is chosen to have the same dimension as Q and if we introduce an (unnormalized) maximally entangled vector, 

I*) = E \9n) ® Ien) , (2.4) 

where {Iffn}} and {|e„)} are orthonormal bases on R and Q, then we can define one-to-one, onto maps of operators 
on Q to vectors on RQ and of superoperators on Q to operators on RQ via 

\Or,q) = Ir®Oq\^) <^=^ {9n,ern\0RQ) ^ (e„JOQ|e„) = (|e„)(e„||OQ) , (2.5) 

Frq =Ir(8)J'q(\^){'^\) <^=4> {gn,em\FRQ\gn',em') = (em|-^Q(|e„)(e„'|)|e™') = (|e„i)(e„|| J"Q||em')(e„'|) . 

(2.6) 

The isomorphism between superoperators on Q and operators on RQ is called the Choi-Jamiolkowski isomorphism [3, 
la]. The second form of Eq. (j2.6|) shows that the left-right positivity of J^q is equivalent to the positivity of Frq. 
Hence, we have the standard result that J-'q is completely positive if and only if it is left-right positive. 

Given any orthonormal basis of operators, {fa}, acting on T-Lq, we can define the unit superoperator acting in the 
left-right sense: 

Iq^E^"0^^=EI^")(^"I- (2-7) 

ot a 

By considering an outer-product basis, { fjk — \ej){ek\ } for some orthonormal basis, { |ej) }, it is easy to check that 

lQ{0)^J2^^'^O^k^^<0)lQ. (2.8) 

The setting of our theorem is that [/ is a joint unitary operator on A and B and as is an operator on B. We 
define a superoperator on A by its action (|l.ll) on density operators. We want £ to map density operators to density 
operators, so it must at least be Hermiticity- and trace-preserving. Requiring that £{pa) be Hermitian does not imply 
that ctb is Hermitian, but it does imply that £{pa) — trs ([^/^^'^'^(o'B-l-cr^)!/^). Since the Hermitian ancilla operator 
2i'^B + crjj) gives the same map as ctb, we can assume that (Tb is Hermitian without loss of generality. Requiring 
that £ be trace preserving implies that as has unit trace. 

We will need the Schmidt decomposition of C/ [3] , 

Ru 
U = J2^n®Bn, (2.9) 

n=l 

where Rjj is the Schmidt rank of U. The Schmidt operators An are nonzero and orthogonal, 

(^„|A„) = tVAiAiArn) = iAn\An)Snm , (2.10) 



and likewise for the Schmidt operators i3„. Sometimes we want to deal with normalized operators, 

i„= , ^" = and B^= , ^" = . (2.11) 

\/(A„|A„) v/(S"l^n) 

We can normalize either set, {A„} or {-B„}, by absorbing the normalization constant into the other set. That U is 
unitary means that 

Ru 



I^U^U^ J2 AiA^(SBiB^, (2.12) 



n,m— 1 

which, by taking the partial trace on A, implies that 

Ru 



/B = ;^^(A„|A„)i3tij„. (2.13) 



dA 

n—l 



Taking the trace on B then gives 



dAdB = Y. (^n|AO(S„|B„) . (2.14) 



We will need the span of the Schmidt operators Bn, 

S = Span{B„ I n= 1,2,--- ,i?[/} , (2.15) 

and the operator orthocomplement of B, 

Ob^{Ob I trB(B^OB) = 0yBeB} = {OB \ trB(t/^Os) = }. (2.16) 

Notice that Ru + dim(C'B) = d|. 

It will also be useful to introduce a set of positive operators associated with B, 

C = {B^B I BeS} . (2.17) 

The set C is not a subspace. It is a cone — i.e., any nonnegative multiple of an element of C is also an element — which 
is a subset of the cone of positive operators on B. It is clearly a closed set. We let Ci denote the set of unit-trace 
elements in C; Ci is a closed and bounded subset of the closed and bounded set of density operators on B. Neither C 
nor Ci is necessarily convex. 

For the remainder of this section, we deal only with system _B, so wc generally omit the subscript B on vectors, 
operators, and superoperators. The most important set in our considerations is the set of vectors \<f)) € 'Hb such that 
\(t>){(j)\ ^ C. Any nonzero multiple of such a \(j)) is also in this set, but we emphasize that the set is not a subspace. For 
this reason, it is most convenient to represent the vectors in this set by their normalized counterparts (state vectors), 
so we define the following set: 

Sb^{ state vectors |0) e Ub \ \'f)){(i>\ iC] . (2.18) 

The set Sb is a property of the Schmidt operators -B„ and, hence, a property of U . 
Since we have restricted Sb to normalized vectors, we can replace the definition by 

Sb^ { state vectors |0) e Ub \ \'f)){4>\ iCi] . (2.19) 

We can think of Sb as the set of pure states not in Ci. Because Ci is a closed set, we can say a bit more. If the 
one-dimensional projector |(/')(0| is not in Ci, then it lies some minimal distance away from all elements of Ci, which 
means that there exists a (5 > such that {(j)\C\(j)) <\ — 5 for all C G Ci. Slightly more generally, we can say that 

state vector |(/)) e Hb such that \(f)){4>\iCi =» 3 a (5 > such that (0|C|(/') < tr(C)(l -5) \/ C eC. (2.20) 

It will turn out to be useful to have available equivalent definitions of Sb- Suppose IV') (01 = B Cz B; then 
{tp\tp)\<j>){(l)\ = B^B = C, which if \ip) is not the zero vector, implies that |0)(V'| <E C. Moreover, if \<j>){(l)\ — C — B^B 



for some C G C (and B G B), the polar-decomposition theorem asserts that B is rank one, and thus B = \ip){(j)\ for 
some normahzed |'0) G 'Hb- Thus we have the fohowing equivalence, 

|V')(0| e B for some nonzero |i/;) e ^3 <^=^ |0)((/)| G C , (2.21) 

and its contrapositive, 

V nonzero \ip) G Hb, H){4'\ i^ <=^ {(/)) {(f>\ (^ C . (2.22) 

Consider now the subspace of Hb generated by the operators in the orthocomplement Ob acting on a particular 
vector \(j)) G Hb'- 

Ob\^) = {0\<j>) I OgOb} . (2.23) 

It is easy to see that Ob\<P) is a subspace; it is the span of vectors On|0), where the set {0„} consists of any operators 
that span Ob- Thus we clearly have diin{0 b\(J)}) < dim(C'B). We denote the orthocomplement of Oslt/)) by {Ob\(J)))j_- 
Because B and Ob are orthogonal subspaces, we have 

|^)(0| e S ^^ V O G Ob, = tr(Ot|i^)((/.|) = ((/.|Ot|^) ^=^ |^) g (Ob|</,))^ (2.24) 

and the contrapositive, 

Wicj^liB ^=^ \^}i{OBm±. (2.25) 

This gives us the following equivalence: 

V nonzero \^Ij) G Hb, W{<f>\ iB <^=^ ObH) =Hb ■ (2.26) 

Thus because of Eqs. (|2.22p and (|2.26p . the following two ways of defining Sb are equivalent to Eqs. p.lSp and 

Sb = { state vectors |0) G Hb \ V nonzero |^) G Hb, \'^){4'\ ^ ^] 

= {state vectors \cj)) G Hb \ Ob\(1>) =Hb} ■ (2.27) 

It is useful to characterize Sb in yet another way. For this purpose, consider the orthonormal Schmidt operators, 
{Bn, n = 1, . . . ,Ru}, of Eq. (j2.1ip . Complete a basis of orthonormal operators by adding a set of orthonormal 
operators, { 0„ | n = Ru + 1, . . . , d^ }, from the orthocomplement Ob- We have 

B.U d% 

I = ^S„0Bt+ J2 OnQOl, (2.28) 

n— 1 n—Ru+1 

and for state vectors \(j)) G Hb, 

Ru <^B 

/ = l(|</))(0|)=^B„|(/.)(0|St+ Y. OM{<t>\Oi. (2.29) 

n=l n=B,u+l 

A state vector |0) is in Sb if and only if the second sum on the rightmost side of Eq. (|2.29p is a full-rank operator, 
i.e., has no zero eigenvalues. This happens if and only if the first sum on the rightmost side of Eq. (|2.29p is strictly 
less than /; i.e., it has largest eigenvalue strictly less than unity. Thus we can characterize Sb in the following way: 



Sb ~ { state vectors \(b) eHi 



1 > 



Rtj 



E^«i<^)('^i^" 



(2.30) 



Here |lM|loo is the operator norm that denotes the largest eigenvalue of vAftM, i.e., the largest singular value of M 
(for positive operators, as here, the largest eigenvalue of M). 



III. MAIN THEOREM 

The content of our theorem is that any rank-one projector not in Ci can be used to construct an unphysical as 
that gives a completely positive map £. Formally, we can state the theorem in the following way. 

Main theorem. Let [/ be a joint unitary operator on A and B and as be a unit-trace Hermitian operator on 
B. Define a Hermicity- and trace-preserving superoperator f on A by £{pa) — trs^t/p^ (8) itbC/^). Ci contains all 
rank-one projectors, i.e., Sb is the empty set, if and only if requiring that f be a quantum operation implies that as 
is a density operator. In symbols, 

V state vector |0b) e Hb, |0b)((/'s| e Ci ■^=^ S's = 

£ a quantum operation => (Tb is a density operator) . (3.1) 



We use P to denote the property of U that Sb is the empty set. ^P is thus that Sb is not the empty set, i.e., that 
there are pure states not in Ci. Any \4>b) G Sb can be used as an eigenstate of ctb with negative eigenvalue in an 
ancilla model that gives a completely positive £. 

The following properties imply P and thus are sufficient to ensure that ctb must be a density operator: 
pi. ( tiBiU^OB) = =^ Ob = 1 , which is equivalent to dim(C'B) = and Ru — dg. 

p2. Ru > d^— da, which is equivalent to dim (Ob) < (Ib- This property implies P because dim (Ob |(/)b)) < dini(OB). 

The following properties imply ^P and thus are sufficient to ensure that there is a nonpositive aB that leads to a 
completely positive £: 

ql. 3 a non-vanishing Hermitian operator Eg such that — trB{Bn'^BBl^) = trBiB^Bn'^B) V n,m. This is 
equivalent to trsiB'SBB'') = for all B G B and thus to trB(CSB) = for all C G C. If a non-Hermitian 
operator satisfies this property, so does its conjugate, thus allowing us to construct a Hermitian operator that 
satisfies the property. That Ss is necessarily traceless follows from Eq. (|2.13l) . Suppose we add /iSs, where /i is 
any real number, to a density operator aB that is diagonal in the same basis as S^, i.e., (t'b = <^b + M^s- Then 
a'g has unit trace, and the matrix G of Eq. (|3.3p — and, hence, the quantum operation £ — remains the same, 
but by making |^| large enough, we can give ct^ negative eigenvalues. All the eigenvectors of S^ are elements 
of^B. 

q2. dim! Spanj B^-^Bn m,n — 1, . . . , Rjj } j < dg. This is equivalent to ql. 

q3. Ru < dB, which is equivalent to dim(OB) > dg — dB- This property implies q2, since dim! Spanj B^Bn } 



< 



Now we are ready to give the proof. 

Proof. We first reformulate the condition for complete positivity of £ in terms of the Schmidt operators on B. 
Plugging Eq. ^^ into Eq. dTl]), we find 

Ru Ru 

£ = / ^ Gnm-An Q A^ = / ^ Gnm\An)[Am\ , (3-2) 

n,m— 1 n,m— 1 

where 

Gnrn = {An\£\A,n) = tr^ (B„(TbS1„) • (3.3) 

Complete positivity of £ is the statement that £ is positive in its left-right action. This is equivalent to the positivity 
of the matrix G, which is that J2n m CnGnmC*m — for all coefficients c„, and this becomes the property 

tTBiBaBB"!) > y B eB. (3.4) 

Wc have reduced our considerations to quantities defined on system B, so we often drop the subscript B for the rest 
of the proof. The complete positivity property of £ can now be written in the compact form 

tr(Ca) > y C eC . (3.5) 



If P holds, C contains all one-dimensional projectors, so the complete positivity condition (I3.5P immediately implies 
that a > 0. We prove the converse in the contrapositive; i.e., we show that if ^P is true, we can construct a nonpositive 
(unphysical) a that is consistent with Eq. (|3.5|) . Since P is not true, there is a state vector \(j)) such that \(f>){4>\ ^ C 
and, hence, by Eq. (|2.20p . a S > such that (0|C|0) < tr(C)(l — S) for all C € C. We lose nothing by also requiring 
that 6 < 1/2. For any such state vector, we can construct a nonpositive a that satisfies Eq. (13.51) by letting \<j)) be an 
eigenvector with a negative eigenvalue — e: 

a = -e|^)(</.| + -l±^(/-|0)((/)|), e>0. (3.6) 

For all C e C, we have 

triCa) = -i±i^ (0|C|^) + -liL tr(C) > (s]±^ - .) tr(C) . (3.7) 

By choosing e < S/{dB — 1 — Sdn), we can ensure that the complete positivity condition p.Sp is satisfied for this 
nonpositive a, thereby completing the proof. ■ 

IV. EXAMPLES 

Example 1. U — Ua ® Ub- In this case, any unit-trace operator gb gives the same (unitary) quantum operation. 
There is a single Schmidt operator, Ub, on B {Ru = 1). The subspace B is the one-dimensional subspace of all 
multiples of Ub, C is the cone of nonnegative multiples of Ib, Ci = { Ib/Ab }, Ob is the (d^ — l)-dimensional space 
of operators of the form UbO, where O is any traceless operator, and Sb contains all state vectors. 

Example 2. The SWAP operator for two systems of the same Hilbert-space dimension: 

u - SWAP - J2 1^^'/^-) ^ (e,'M = J2 le'»)(e.l ® \fj){M ■ (4.1) 

Here the vectors \ej) and \fk) make up any orthonormal bases for A and B. The SWAP operator acts on product 
states as UpA (^ ctbU^ = (ja® Pb, so it is clear that ctb must be a density operator. Equation (|4.ip is a Schmidt 
decomposition of SWAP. Schmidt operators on B can be taken to be the outer products \fj){fk\- These operators 
span the operator space on _B, so B contains all operators, C is the cone of positive operators, Ci is the set of density 
operators, Ob consists only of the zero operator, and Sb is empty. 

Example 3. Two qubits with a joint unitary 

C/ = e''^®^^/2 = /®/cos(6'/2)+iX(g)Xsin(6'/2) , < 6* < tt . (4.2) 

This joint unitary has Schmidt rank 2. For the unitary of Eq. (14. 2p . the two Schmidt operators for B can be taken to 
be / and X, Ob = Span{y, Z}, C consists of all operators of the form r{I + sX), where r > and \s\ < 1, Ci is the 
subset of C with r = 1/2 (i.e., the density operators whose eigenstates are the eigenstates of X), and Sb contains all 
state vectors except the eigenstates of X. 

If we write as as ctb = i(/ + s • cr) = ^{I + s^Xb + SyYB + s^Zb), where cr is the vector of Pauli operators on B, 
it is straightforward to show that the superoperator (jl.ll) on A is 

£ = cos\9/2)I Ql + is^ sin(6i/2) cos(6'/2)(X Q I ~ I Q X) + sm\e/2)X QX = p+V QV'' + p-XV V'' X , (4.3) 



where p± — ^^(1 ± T), T — J cos^O + s^ sin^0, V = e*'^"/^, and tana — s^ tan 6*. This map is a convex combination 

of two orthogonal unitaries, V and XV , and is completely positive if and only if T < 1, i.e., \sx\ < 1. The point here 
is that Sy and s^ can be anything without changing £ . Complete positivity only constrains s^, in agreement with our 
conclusion that the only state vectors not in Sb are the eigenstates of Xb ■ 

Any two-qubit joint unitary has Schmidt rank 1, 2, or 4 8]. A Schmidt-rank- 1 unitary is a product unitary and 
thus is covered by Example 1. A Schmidt-rank- 2 unitary is equivalent under local unitaries to Eq. (|4.2p and thus 
covered by Example 3. A Schmidt-rank-4 unitary has full Schmidt rank; Ci is the set of all density operators, and 
Sb is the empty set. Since Ci is convex in all three cases, we can conclude that for two-qubit unitaries, C and Ci are 
always convex sets. 

Example 4- A qutrit and a qubit with joint unitary 

[/= |0)(0|«)/+|l>(l|«)X + |2)(2|(g)Z. (4.4) 



This joint unitary has Schmidt rank 3. The Schmidt operators on B can be taken to be I, X, and Z, Ob consists of 
ah multiples of Y, C is the cone of positive operators, Ci is the set of all density operators, and Sb is empty. That U 
satisfies p2 makes clear that Sb is empty. 

For the remaining examples, we refer to an orthonormal basis, { |j), j = 0, . . . ds — 1 }, on system B, we denote 
the rank-one outer products by Tjk = \j){k\, and we use the following operators on B, defined for j < k: 

Ijk = Tjj + Tkk , (4.5) 

^jk = Tjj - Tkk , (4.6) 

^jk = Tjk + Tkj , (4.7) 

Yjk = ~i{Tjk - Tkj) . (4.8) 

These are the Pauli operators for the two-dimensional subspace spanned by \j) and |fc). The following operators make 
up a set of dg linearly independent operators, which span the space of operators on B: the unit operator /, and the 
operators Zqj, j = 1, . . . , ds — 1, span the space of diagonal operators; the operators Xjk, numbering dB{dB — l)/2, 
and the operators Yjk, numbering dB{dB^ l)/2j span the space of off-diagonal operators. These operators are pairwise 
orthogonal, except that trlZQjZok) = 1 + Sjk- 

We can introduce an associated set of dj — 1 unitary operators: 

,dB-l, (4.9) 

(4.10) 
(4.11) 

The unitary Uj is diagonal, with ones everywhere on the diagonal, except in the jth position, which has a —1. The 
operators Vjk and Wjk act like the Pauli X and Y in the two-dimensional subspace spanned by \j) and \k), and act 
like the identity operator on the orthocomplement of this subspace. Together with the identity operator /, these 
operators are linearly independent and span the space of operators on B, but they are not orthogonal. To see that 
they span the space, one notes that the diagonal operators satisfy 



Uj^I- 


- loj + Zqj , 


.? = !,• 


^,k^I- 


- Ijk + Xjk , 


J < k, 


'^,k^I- 


- Ijk + Y,k , 


j < k. 



^i, = ^i'-Vi). j = l.^^-,dB-i. (412) 



These expressions allow us to write any diagonal operator in terms of / and the unitaries Uj . Writing the diagonal 
operators I — Ijk in terms of/ and the unitaries Uj, we find that the operators Xjk and Yjk, which span the off-diagonal 
subspace, can be expanded in terms of our set of unitaries. 

Example 5. This example shows that P can be satisfied even though Ru — Sds — 2 is linear in dB- The idea 
is to introduce a unitary U such that the set B contains the operator subspace spanned by the operators tqj for 
j = 0, . . . , (is — 1, i.e., contains all operators of the form 

B^XJ2 b*roj , (4.14) 

where A > and the expansion coefficients bj arc normalized to unity. This means that C contains all operators of 
the form 

B^B = X^J2^,blT,k = A2|0)(</)| , (4.15) 

j,k 

where 

m^Y.^j\j) (4.16) 

is a state vector in Hb- Thus Ci contains all pure states, and Sb — 0- 
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To do this, we assume that system A has dimension Sds — 2, with an orthonormal basis of vectors |0), \rj), \sj), 
and \tj), where j = 1, . . . , d^ — 1, and we introduce the foUowing controlled unitary on A and B: 

t/ = |0) (0| / + ^ {\T,){r, \®U, + \s,) {s, I Voj + \t,){t, I Wo,) . (4.17) 

i=i 

Because the projectors corresponding to the orthonormal basis on A arc orthogonal and, hence, linearly independent, 
B is the span of the target unitaries, and the operator subspacc Ob consists of all operators that are orthogonal to 
all of the target unitaries. By the same argument as above, we then have 

B = Span{ Too; Tjj, tqj, and tjq, for j = 1, . . . , ds - 1 } . (4-18) 

As noted above, in this situation, Ci contains all pure states and thus P is satisfied. 

We can go a bit further and characterize the entire set Ci . An arbitrary element of B has the form 

B = ^J2 ^j^oj + J2 ^'M^J^ + '^*no) > (4-19) 

where A and /ij are real and nonnegative, the expansion coefficients hj are normalized to unity, and jc^p + \dj\'^ — 1. 
The general element of C looks like 

C = A2|^)(0| + ^ ^i]{c,\Q)+d,\J)){c*{Q\+d*{3\) , (4.20) 

i=i 

where |(/)) is the state vector of Eq. (|4.16p . Thus the general element of Ci is a convex combination of an arbitrary 
pure state \(j)){(t>\ and the pure states in the sum in Eq. (|4.20p . 

We can use this example to show that Ci (and C) are not generally convex sets. If Ci were a convex set, then since 
it contains all pure states, it would have to contain all density operators, but it is easy to find density operators not 
in Ci for ds > 4. For this purpose, we restrict attention to operators that are orthogonal to |0). The general form of 
an element of Ci of this sort is 

C = \' £M:b-)(fc|+EA^'|j)01- (4.21) 

Let n be the projector onto the subspace orthogonal to |0), and let S be any traceless, Hermitian operator that is 
orthogonal to |0), scaled so that its most negative eigenvalue is —1. Then a — {Il + e'S)/{dB — 1) is a density operator 
for < e < 1; when e = 1, ct is on the boundary of the space of density operators. For this density operator to be in 
Ci, we must have X^bjb^. = eY.jk/{dB — 1). When ds > 4, we can let E12 = S23 = — S13 be real and positive; then 
there is no way to choose the phases of 61, 62, and 63 to satisfy the condition if e 7^ 0, thus showing that a is not in 
Ci for 0< e< 1. 

It is worth spelling this out in three dimensions (ds = 3), where the target unitaries on B have the following matrix 
representations: 

/l 

/=|oio|, ;7i=(o-io|, ;72=oi o 

\0 -1, 







V^oi = 1 , ^02 = 1 , Woi^ \i 0\ , W02 = 1 . (4.22 



We can see explicitly in this case that B has the form given in Eq. (I4.18P . The orthocomplement Ob , being orthogonal 
to all the target unitaries, is 



0\ 


/GOO 


1, 


-i 


1 0/ 


\o i 



Ob = Span M001,00-iV= Span MOOlkOOOV. (4.23) 



0\ 


/O 


1, 





0/ 


\o 1 



Any operator Ob € Ob has all zeros in the its first row, so {0\Ob\(I>) = for any vector |(/)), implying that Sb — ^■ 
This way of showing that Sb is the empty set, by showing that all elements of Ob have all zeroes in the first row, 
works for any value of cIb ■ 

Example 6. It is easy now to formulate an example where P holds and Rjj — 2dB- To do so, we consider the controlled 
unitary 

ds-l 

U^\0){0\®I+\n){n\®U[+ J] (|.s,)(s,| ® Fo, + |i,)(t,| ® P^^,) . (4.24) 

where the unitaries on B are defined by 

U[=I-Ioi- Zoi , (4.25) 

Vo, =1- hj + Xo, , (4.26) 

Wl,^^I-h,-iY^,. (4.27) 

It is easy to convince oneself that the target unitaries are linearly independent and thus that Ru — 2dB- Moreover, 
since tqo = (/— U[)/2 and tqj — (Vqj — Wqj)/2 = (Xqj +iloj)/2 for j = 1, . . . , ds — 1, we are in the situation discussed 
in Example 5, so Sb = 0- 

Example 1. In this example, we demonstrate that ^P is possible when Ru — dg ~ cIb, thus showing that condition p2 
is optimal. We assume that system A has dimension d\ — dB, and we use the controlled unitary 

U^\0){0\®I+Y. I^j) (^J- 1 ® f^J + £ £ (l^Jfe) (*Jfe I ® ^1^ + l*Jfe) ^^ok I ® M(,-fe) + £ l^oj) (toj I ® Waj , 

(4.28) 

where the indicated vectors on A make up an orthonormal basis. As in Example 5, B is the span of the target unitaries 
and thus can be written as 

B = Span{ Too; T^k, for j, fc = 1, . . . , ds - 1; Foj, for j = 2, . . . , ds - 1 } , (4.29) 

and the orthocomplement of B is 

Ob = Span { roi, no; Xqj, for j = 2, . . . , ds - 1 } . (4.30) 

Letting |0) = |0) + |1), we clearly have Ob\4') — "^s and thus -iP. Indeed, it is easy to see that the state vectors not 
in Sb are those that are orthogonal to |0) or to |1). 

To introduce our last example, we consider the Weyl (or generalized Pauli) unitary operators 

Bm^X'^Z', k,l = 0,...,dB~l, (4.31) 

where X and Z are the shift and phase operators in the standard basis, i.e., X\j) — \{j + 1) mod ds) and Z\j) ~ 
e'^'^-'/'^^ |j). The Weyl operators (|4.31|) are pairwise orthogonal. The eigenstates of X are the Fourier-transformed 
basis states 



, dfl— 1 

liy = -^= y e'2-jfe/dB |j-^ (4 32) 



i=o 



We have X\k) = er''^^^l'^^\k) and Z\k) = \(k + 1) mod rfs). 

We assume the existence of a fiducial state vector |0) such that the state vectors \4>ki) — Bki\<t>) form a SIC-POVM, 
i.e., 

\{M<l>k'i')\^ = :r^ when kl^k'l'. (4.33) 

The existence of such a fiducial vector for all ds was conjectured by Zauner [9|; the conjecture has been demonstrated 
for ds = 1, • • ■ , 15, 19, 24, 35, and 48 and numerically up to at least dg = 67 (see |T0[ and references therein). Notice 
that any vector in the SIC-POVM could be used as the fiducial vector. 



Consider now the class of controlled unitaries with Schmidt decomposition 

Ru 



C/==^|n)(n|«B„, 
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(4.34) 



n=l 



where the operators -B„, indexed by the single index n, comprise a subset of Ru elements drawn from the Weyl 
operators (|4.3ip . The form (|4.34p is a Schmidt decomposition of U, and the Weyl operators i?„ are Schmidt operators 
on system B. 

According to (|2.30p . the fiducial state vector \(j)) is in Sb if and only if 



Here we define 



Ru Ru 

5^B„|(/.)(0|Bt = ^ !</.„) (0„ 



Ru 
n=l 



l$lloo < dp 



(4.35) 



(4.36) 



to be the sum of the projectors onto the subset of SIC-POVM vectors. Our task is to bound the operator norm of <I> 
in terms of Rjj. 

To do so, we consider a purification of $ into an i?(7-dimensional system S: 



\<\>sb) ^ £ 



(4.37) 



n=l 



Since |ct)sB) is rank one, the spectrum of $ is the same as the spectrum of the reduced matrix of |4'sb)(4^sb| on S, 
which is the Gram matrix of the vectors loin): 



Ru 



trs(|4)ss)(4)5s|) = ^ {<f>rn\<l>n)\n){m\ 



(4.38) 



n,m— 1 

Applying the condition (|4.33l) that the vectors |(/)„) form a SIC POVM gives 

Ru 



trB(|4)5B)(4>5B|) = ^|n)(n| + -=i== J2 



\n){m\ 



n,r?i— l....,i?L 



1 \ ^ Ru 

1- , ]Is+ , Y e^^— |n)(m 



(4.39) 



(4.40) 



n,rn— 1 



with 9„m = arg(((/)„|0m)) = ~Omn- Since the first term in Eq. (|4.40p is proportional to the identity on S, we only need 
to bound the norm of the second term, which we do by considering an arbitrary normalized vector |c) = X]n=i ^^ \n): 



Ru 

E 



|n)(7TT- 



^n,7n— 1 

Putting everything together, we obtain 



=. g e''-c*„cm< £ |c„||c™|= (g|c„|) <Ru. 



||$||oo = ||trB(|4)Si3)(4>SB|)L<l + 

The right-hand-side of Eq. (|4.42p is strictly smaller than d^ whenever 



Ru-l 
\/dB + 1 



Ru < iVdB + HdB - 1)1 
Here \x\ denotes the smallest integer larger than or equal to x. 



(4.41) 



(4.42) 



(4.43) 
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We conclude from the condition (I4.35|) that Sb ^^ for Ru < \\/dB + 1(^5 — 1)1, regardless of the specific choice 
of the Rjj Weyl operators i3„ in Eq. (j4.34|) . In particular, Ss contains the fiducial vector |0). We now examine in 
detail a particular case of the joint unitary (|4.34|) . 

Example 8. Consider the controlled unitary (|4.34l) with the specific choice of Ru = 2dB — 1 Weyl operators 



{/, Z, . . . , Z''^ 



X. 



^^dB-1} 



(4.44) 



For dB > 4, condition (|4.43p is satisfied; hence the fiducial state vector is contained in Sb- Moreover, the special 
Fourier structure of this case makes it easy to determine all the elements of the set Sb for rfs > 3 (the case d^ = 2 is 
equivalent to our Example 4). To do this, we find the set of state vectors contained in Ci; Sb is the complementary 
set. For dB > 3, we now show that the state vectors in Ci are the eigenstates \j) of X and the eigenstates \j) of Z. 
The general element of B, 






(4.45) 



has the matrix representation in the standard basis. 



B = 



/ Co "ndB-i VdB-2 

m Ci VdB-i 
m vi C2 

\VdB-l VdB-2 VdB-S. 



Cdfl-i/ 



(4.46) 



where 



C.=A+^ 



^i2TTJk/ di 



Ck 



(4.47) 



k=l 



The matrix B is rank one if and only if all the order-two minors (determinants of 2 x 2 submatrices) are zero. 

We first note that minor condition implies that if any of the rjjS is nonzero, then all the matrix elements of B — 
hence, all the rjjS — are nonzero. We deal with this case first. The condition on minors can then be restated as the 
requirement that the ratio of matrix elements in the same row of different columns is independent of row, and these 
reduce to the following conditions: 



_rn_ 
Cj+i 



m 
m 






VdB-l 

VdB-2 



C, 



VdE 



j = 0,...,dB- 1- 



(4.48) 



For dB = 2, this is a single condition, r/f — CoCi = ^'^ ^ Cij as mentioned above, this takes us back to Example 4, 
where Sb — 0- For dB > 3, we have immediately that all the C s are equal and thus that Q — ior j — 1, . . . , dB — I 
and C,- = A for j = 0, . . . , (is — 1. Now we can find the solutions of conditions (|4.48p . labeled by fc = 0, . . . , d^ — 1: 



Vj 



\^^i2TTJk/dE 



(4.49) 



This gives Bk — dB^\k) {k\. We now turn to the case where all of the rjjS are zero. In this situation, the minor 
condition immediately tells us that only one diagonal element of B can be nonzero. We can also deal efficiently 
with this case by moving to the Fourier basis, which switches the roles of X and Z and of the rjjS and CjS; then one 
similarly concludes that aU the ry^s are zero and that Cj = \Q-^2-Kjk/dB g^jj^ f^^ _ dsXSjk, fc = 0, . . . , d^ — 1, giving 
Bk — dB^ \k){k\. This shows that the pure states in Ci are the eigenstates of X and Z. All other state vectors are in 
Sb- 



INDIRECT TOMOGRAPHY 



Suppose that Sb is not the empty set so that, as we have learned, there are unphysical ancilla "states" ctb ^ 
such that the map (jl.ip is a proper quantum operation. One can ask whether this happens because the unphysicality 
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of the ancilla is "hidden" and irrelevant. In the example given in the Introduction, i.e., a joint unitary U = Ua ® Ib, 
the ancilla does not interact with the system; the evolution the system undergoes is independent of the "state" of the 
ancilla. We will now show that this is not always the case: a counterintuitive situation can arise, where (i) the joint 
unitary U can lead to an operation on the system even in cases where the "state" of the ancilla is unphysical, i.e., 
Sb 7^ 0, yet (ii) at the same time, U allows the tomographic reconstruction of the "state" of the ancilla. 
By this we mean the following. Recall Eqs. (|1.1|) . p.2p . and p.3|) . 

Ru 
£a{pA) = tTsiUpA ® (TbU'^) = Y^ GnmAnPAAl, , (5.1) 

n,rM— 1 

where Gnm = {An\£a\Am) = iiB{Bn<^BBJ^), and the operators An and B„ are Schmidt operators for U, as in Eq. (12.91) . 
We can use standard process tomography to determine £^ from measurements of an informationally complete POVM 
on the output states obtained from an informationally complete set of input states. We say that U allows indirect 
tomography on B if this knowledge of £cr uniquely determines cfb ■ This is equivalent to saying that the matrix Gnm 
determines as- 

To allow such indirect tomography on the ancilla state, all we need is that different operators as map to different 
matrices Gnm- This map being linear, we require that trB{BnOBB^^) — 0, ioi n,m — 1, . . . Rjj, implies Ob — or, 
equivalently, that the operators B^^Bn span the space of operators on Hb- Thus we reach our theorem on indirect 
tomography. 

Theorem. Indirect tomography. The unitary U with Schmidt decomposition U = X^n^i An ^ Bn allows indirect 
tomography if and only if 

dim! Span{ BJ^Bn \ m,n = 1, . . . , Ru } j — dg . (5-2) 

We want to exhibit unitaries that allow indirect tomography even if Sb y^ 0- We do not have to go far: Examples 7 
and 8 of Sec. IIVI have exactly this property. In Example 8, it is not hard to see that the set {BJ^ Bn \ m,n = 
1, . . . , 2dB — 1 } comprises all Weyl operators (up to irrelevant phases), which means that the unitary (I4.34p allows 
indirect tomography. 

By definition, when an interaction unitary U allows indirect tomography there is only one "state" ctb of the ancilla 
compatible with the evolution f^- Thus the situation incarnated by Examples 7 and 8 clarifies that the complete 
positivity of the evolution £^ undergone by A, in spite of the lack of physicality of ctb, cannot be interpreted as saying 
that there is some other physical state a'^ > of the ancilla such that So- — £a' ■ 

VI. CONCLUSION 

We have introduced, answered, and explored the question of what properties a joint unitary operator, acting on 
a system and an ancilla, must have in order that a quantum operation on the system, based on this joint unitary, 
requires that the ancilla state be physical. The answer uncovers structures, not hitherto appreciated, on the space of 
states of the ancilla, structures clearly connected to the question of when a superoperator is a completely positive. 
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